Percolation quantum phase transitions in diluted magnets 
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We show that the interplay of geometric criticality and quantum fluctuations leads to a novel 
universality class for the percolation quantum phase transition in diluted magnets. All critical 
exponents involving dynamical correlations are different from the classical percolation values, but in 
two dimensions they can nonetheless be determined exactly. We develop a complete scaling theory 
of this transition, and we relate it to recent experiments in La2Cui_ p (Zn,Mg) p 04. Our results are 
also relevant for disordered interacting boson systems. 



Geometric criticality and quantum criticality are two 
distinct phenomena leading to universal scale invariant 
correlations. Both get combined in randomly diluted 
quantum magnets. Site or bond dilution defines a per- 
colation problem for the lattice with a geometric phase 
transition at the percolation threshold . Quantum fluc- 
tuations of the spins coexist with geometric fluctuations 
due to percolation. In this Letter, we address the ques- 
tion of whether or not these quantum fluctuations fun- 
damentally change the percolation phase transition. 

The generic phase diagram of a diluted magnet as func- 
tion of impurity concentration p, quantum fluctuations 
strength g, and temperature T is shown in Fig. ^ In 
the absence of both thermal and quantum fluctuations 
(g = T = 0), magnetic long-range order survives for all 
concentrations p < p c , where p c is the classical (geomet- 
ric) percolation threshold. If p > p c , no long-range order 
is possible as the system is decomposed into independent 
clusters. For p < p c , magnetic order survives up to a 
nonzero temperature T c (p); but at p c , it is destroyed im- 
mediately by thermal fluctuations, T c (p c ) = 0,0,0,0- 

Zero temperature quantum fluctuations are less effec- 
tive in destroying magnetic order. Even at p c , the perco- 




FIG. f : Schematic phase diagram of a diluted magnet with 
impurity concentration p, temperature T, and quantum fluc- 
tuation strength g |2J. There is a quantum multicritical point 
(big dot) at (p c ,g*)- The zero temperature quantum phase 
transition across the dashed line is the topic of this Letter. 



lating cluster can remain ordered up to a nonzero fluctu- 
ation strength. This was found in the diluted transverse 
field Ising model 0, Q where a multi-critical point with 
transverse field strength g* and p — p c emerges. At this 
point, both quantum and geometric fluctuations diverge. 
In O(N) systems (N > 2), early results suggested that 
quantum fluctuations always destroy long-range order at 
Pc 03- E El- However, Sandvik ^1 showed that the 
percolating cluster in a 2d diluted Heisenberg quantum 
antiferromagnet (QAFM) is ordered at p c (see also Ref. 

If quantum fluctuations are enhanced, e.g., by an 
interlayer coupling in a bilayer QAFM Il5| , a multi- 
critical point at (p c ,g*) arises, similar to the transverse 
field Ising case. The behavior at the generic quantum 
phase transition occurring for g > g* , p < p c (solid line 
in Fig. ^) and at the multi-critical point at (p CJ ff*) has 
been investigated in the past |l4Lll5lll(Tl |. In contrast, the 
percolation quantum phase transition at p c and g < g* 
which can be observed, e.g^ in diluted 2d QAFMs such 
as La2Cui_p(Zn,Mg) p 04 |ljfl has received less attention. 

In this Letter, we show that the interplay of geometric 
criticality and quantum fluctuations leads to a novel uni- 
versality class for this percolation quantum phase tran- 
sition. Even though the transition is driven entirely by 
geometric criticality of the underlying classical percola- 
tion problem, quantum fluctuations enhance the singu- 
larities in all quantities involving dynamic correlations. 
As a result, the susceptibility exponent 7 as well as a, 
5 and r\ differ from the classical percolation values, and 
the dynamical exponent z is nontrivial. In contrast, the 
order parameter and correlation lengths exponents are 
classical, /3 — C and v = v c [l^- ^ n the remainder of 
this Letter we sketch the derivation of these results via 
a scaling approach (we have also performed detailed ex- 
plicit calculations) and calculate critical exponent values 
in 2d and 3d. In addition, we make predictions for vari- 
ous experimentally relevant observables. 

Our starting point is a general O(N) quantum rotor 
model (N > 2) on a d-dimensional hypercubic lattice 
with nearest neighbor interactions as arises as low-energy 
theory of a Heisenberg quantum antiferromagnet |l9| . 
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FIG. 2: For g < g* , a percolation cluster acts as a compact 
object. The spins are correlated in space but collectively fluc- 
tuate in imaginary time r. 



The action reads 
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Here Si(r) is a iV-component unit vector at lattice site 
i and imaginary time r. The quenched random variable 
€i describes the site dilution. It takes the values and 
1 with probabilities p and (1 — p), respectively. The pa- 
rameter g in the dynamic term measures the strength 
of the quantum fluctuations [2(J, to n are bosonic Mat- 
subara frequencies, and zq is the dynamical exponent of 
the clean system. We are mostly interested in the case 
zq = 1 (undamped rotor dynamics) but our results also 
apply for < z < 2 where z — > corresponds to the 
classical limit while zq = 2 is the case of overdamped spin 
dynamics (occurring in itinerant magnets |2lll22|). 

Let us start our discussion of the percolation quantum 
phase transition at p — p c and g < g* by considering 
order parameter m and correlation length £. Long range 
magnetic order can only develop on the infinite percola- 
tion cluster; all finite clusters are completely decoupled. 
Since for g < g* the infinite cluster remains magnetically 
ordered for all p < p c , the total order parameter is pro- 
portional to the number of sites in the infinite cluster, 
m ~ Poo{p) ~ [p c — pY c ■ Thus, P = p c . To determine 
the magnetic correlation length £ we note that the corre- 
lations cannot extend beyond the connectedness length £ c 
of the percolating lattice because the rotors on different 
clusters are decoupled. On the other hand, for g < g* , 
all rotors on the same cluster are strongly correlated in 
space even though they collectively fluctuate in time (see 
Fig. |2| . We thus conclude £ ~ £ c and v = v c . 

To study quantities involving dynamic correlations, we 
develop a scaling theory of the free energy, contrasting 
classical and quantum cases. Consider a single percola- 
tion cluster of finite size (number of sites) s. In a classical 
magnet for T — > 0, the free energy of such a cluster has 
the scaling form F s c (h, T) = T$ {hs/T) 7 where h is the 
field conjugate to the order parameter and $ (x) a di- 
mensionless function. The susceptibility x c s follows from 



taking the second derivative with respect to h; this gives 
~ s 2 /T as expected. Let us now consider the quantum 
case. For g < g* , a quantum rotor cluster of size s can 
be treated as a compact object that fluctuates in time 
only (see Fig- El) an d is described by a 0+1 dimensional 
nonlinear sigma model (NLSM) with the action [2^, 



A a = sA dyn [S] + hs J drS {1) (t) 



(2) 



Here, is one component of S. Dimensional analysis 
shows that the free energy of a cluster of size s is 



F s (9,h,T) 



g v s~ 



°<$>(hs 1+(p g-' p ,Ts' p g- v ) (3) 



with exponent ip = zq/{2 — zq) (for zq < 2). This re- 
sult can be obtained by rescaling imaginary time, t — » 
gfg-Vf-^ which eliminates the factor (s/g) in the dynamic 
part of the action, and leads to a dimcnsionlcss field 
strength (hs 1+v / g v ) and temperature (Ts v /g v ). The 
resulting order parameter susceptibility Xs ~ s 2+v of a 
quantum spin cluster diverges more strongly with s than 
that of a classical spin cluster. We have obtained the 
same results from an explicit large- N calculation and 
from a renormalization group analysis of the NLSM (j2J 
at its strong coupling fixed point 24, 26] 

After having analyzed a single cluster we turn to the 
full percolation problem. The total free energy can be 
obtained by summing the cluster free energy F s over all 
percolation clusters. From classical percolation theory, 
we know that close to p c the number n s of occupied clus- 
ters of size s per lattice site obeys the scaling form 



(*) 



(4) 



Here t = p — p c , and t c and a c are classical percola- 
tion exponents. The scaling function f(x) is constant 
for small x and drops off rapidly for large x. All clas- 
sical percolation exponents are determined by t c and a c 
including v c = (r c — l)/(da c ) and the fractal dimension 
Df — d/(r c — 1) of the percolating cluster 0. 

We first briefly recapitulate the percolation transition 
of a classical magnet at T = 0. The total free energy 
is given by the sum F c (t,h,T) = J2 S n * (*) F s ( h > T )- 
Rescaling s — ► s /b D f and taking the limit T — ► yields 



F c (t,h) = b- d F< 



hb u ' 



(5) 



The critical behavior of thermodynamic quantities can 
be obtained by taking the appropriate derivatives. All 
critical exponents coincide with the classical (geomet- 
ric) percolation exponents which is expected because in 
a classical magnet at T = geometric fluctuations are 
the only fluctuations present [jj, y, y| ■ 

We next generalize these ideas to the percolation quan- 
tum phase transition in the quantum rotor model 
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for g < g*. Averaging the free energy of a quan- 
tum rotor cluster over all cluster sizes, F (t, h, T) = 
^2 s n s (t) F s (ffj h,T), and rescaling s — > s/b Df yields 

F(t,h,T) = b- ( - d+z '> F (tb l ' v ,hb {D f +z \Tb z ^j (6) 

with the correlation length exponent identical to the clas- 
sical value, v — v c . z — tpDf plays the role of the dy- 
namic critical exponent. We have verified this interpre- 
tation by calculating the correlation time £ r and showing 
£ T ~ £ z . Comparing the scaling form of the free energy 
© with its classical counterpart we notice that the 
quantum theory is obtained from the classical one, as 
usual, by replacing d by d + z. In addition, one should 
replace Df by Df + z. This occurs since the magnetic 
field couples to the spins at all " imaginary time points" . 

The critical exponents at the percolation quantum 
phase transition can be obtained by taking the appro- 
priate derivatives of the free energy. All exponents are 
fully determined by two classical percolation exponents 
(e.g., the fractal dimension Df and the correlation length 
exponent v — v c ) as well as the dynamic exponent z (the 
classical exponents are recovered by setting z = 0): 



2 -a = (d + z)is 

(3 = (d-Df)u 

7 = (2D f -d + z)v 

S = (D f +z)/(d-D f ) 

2-i] = 2D f -d + z . 



(J) 
(8) 
(9) 
(10) 

(11) 



The exponents a, 7, 5, and rj are modified compared to 
classical percolation, while f3 depends only on d— Df and 
is thus unchanged. The exponents fulfill the usual scaling 
laws, and hyperscaling is valid for d < 6. This means, 
the upper critical dimension is identical to its classical 
value, it is not changed by the substitution d —>■ d + z. 

The percolation exponents are known exactly in 2d , 
as is the dynamical exponent z = Dfip = DtZo/(2 — zq). 
Thus, all exponents for the quantum phase transition can 
be determined exactly. For higher dimensions, v and Df 
are known numerically with good accuracy yj. Table I 
shows exponent values for zq = 1 in 2d and 3d. The de- 
pendence of the critical exponents on zq is shown in Fig. 
El In the limit zq — > 0, the quantum exponents become 
identical to the classical percolation exponents, and the 
dynamical exponent z vanishes. For zq — > 2, the dynam- 
ical exponent z as well as 7 and S diverge, indicating 
activated scaling at an infinite-randomness critical point 
as in the Ising case |8|. 

The scaling form JU) of the free energy also determines 
the temperature dependence of the order parameter sus- 
ceptibility x(i,T) = T-tA^)9 x (^7T) and the heat 
capacity C = T(d 2 F/dT 2 ) = T d / Z Q C {t vz /T) with scal- 
ing functions XjC (x). At p — p c (or for p— p c < T 1 /( I/Z )) J 
this yields X (T) oc T^/^ and C (T) cx T d ' z . Spin 







2d 




3d 


Exponent 


classical quantum 


classical quantum 


a 


-2/3 


-115/36 


-0.62 


-2.83 


P 


5/36 


5/36 


0.417 


0.417 


7 


43/18 


59/12 


1.79 


4.02 


5 


91/5 


182/5 


5.38 


10.76 


V 


4/3 


4/3 


0.875 


0.875 


V 


5/24 


-27/16 


-0.06 


-2.59 


z 




91/48 




2.53 



TABLE I: Critical exponents at the classical and quantum 
percolation transition for d — 2 and 3 in case of zo — 1. 
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FIG. 3: Exponents 7, S, and z as functions of zq in 2d. The 
limiting cases zo — > and zo — > 2 are explained in the text. 



wave type excitations within a cluster have typical ener- 
gies to> sw oc s~ z °/ Df that are large compared to the energy 
gap wo oc s~ v = s~ z l Ds of the cluster because z > zq (for 
zq = 1). Thus, spin waves that were shown to be impor- 
tant away from p c [27} , are subleading at p c as can also 
directly be seen from the heat capacity contributions. 

A scaling theory for the dynamic susceptibility can 
be derived following similar arguments. We obtain 
x(t,u,T) = b~f/ u x(tb 1/u \uob z ,Tb z ) which at p c yields on 
the real frequency axis: 



Imx (u 



i(r 



w-T/C^n (w/T) . 



(12) 



with dimensionless scaling function Q (x) . Eq. (|12fl im- 
plies that low-temperature inelastic neutron scattering 
experiments at the location of the Bragg peak should see 
a sharp increase ~ uj~~ < I( uz } in the scattering intensity for 
lu — > 0, while aW^ z )lmx (u> + i0 + ) obeys w/T-scaling. 

In the remaining paragraphs we discuss experiments, 
simulations and the generality of our results. The NLSM 
(JIJ is the exact low-energy theory of a dimer-diluted bi- 
layer Heisenberg QAFM or in d > 2. that 

of an appropriately foonrf-diluted system. One candidate 
material in 3d is (Tl,K)CuCi3 [2||; interesting quasi-2d 
compounds are SrCu2 (1303)2 or BaCuSi206, where suit- 
able dopants remain to be found. In addition, the ac- 
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tion Q describes Josephson junction arrays with random 
missing junctions and diluted bosons in optical lattices. 

If magnetic ions are replaced by nonmagnetic ones 
in a QAFM, e.g., Cu by Zn or Mg in La 2 Cu0 4 or 
YBa2Cu306, impurity induced moments arise which 
leads to random Berry phases in the NLSM [2^. These 
moments produce Curie type contributions that dom- 
inate the low-temperature behavior in the disordered 
phase p > p c . However, at p = p Cl their contributions are 
subleading compared to the critical singularities. Thus, 
our theory also describes the percolation quantum phase 
transition in site-diluted systems, and we propose to 
study QAFMs such as La2Cui_ p (Zn,Mg) p 04. Recent 
measurements 17] indeed gave a nontrivial dynamical ex- 
ponent (with a somewhat smaller value of z ~ 1.4). The 
measurement of several observables may help to identify 
the scaling regime and possibly refine the value of z. An 
issue that may limit the applicability of our theory very 
close to p c is dilution induced frustration. 

We now turn to simulations, focusing on the dynam- 
ical exponent z. Early series expansions gave values 
in the range of 1.5 to 2. Kato et al. 01 obtained spin- 
dependent values between 1.3 and 2.5 but under the as- 
sumption that the critical percolation cluster is not long- 
range ordered. Recent simulations of a bilayer QAFM 
[I3| and of a bond-diluted 2d QAFM :30j gave values of 
about 1.9, in excellent agreement with our theory. These 
two papers also gave arguments for z = Df as obtained 
in our theory for Zq = 1. 

Finally, we discuss the generality of our theory. To de- 
rive the nonclassical scaling theory © of the percolation 
quantum phase transition we only relied on the fact that 
the energy gap luq of a cluster depends on its size s via 
a power law, luq cx s~ v with <p > 0. If tp diverges, which 
happens in our model for zq — 2, the relation between the 
gap and s becomes exponential. The transition is then 
of infinite-randomness type as in the diluted transverse 
field Ising model These observations agree with a 
general classification of dirty phase transitions according 
to the effective dimensionality of the droplets or clus- 
ters |2f| . For zq < 2, the clusters (which are finite in 
space but infinite in imaginary time) are below the lower 
critical dimension d~ of the problem resulting in power- 
law scaling (and exponentially weak Griffiths effects) . In 
contrast, for zq = 2, the clusters are exactly at d~. This 
leads to activated scaling and power-law Griffiths effects. 
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